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Abstract 

The notion of spectral flow has given new insight into the motion of vortices in superfluids and super- 
conductors. For a BCS superconductor the spectrum of low energy vortex core states is largely determined 
by the geometric optics limit of Andreev reflection. We use this to follow the evolution of the states when 
a stationary vortex is immersed in a transport supercurrent. If the core spectrum were continuous, spectral 
flow would convert the momentum flowing into the core via the Magnus effect into unbound quasiparticles 
— thus allowing the vortex to remain stationary without a pinning potential or other sink for the inflowing 
momentum. The discrete nature of the states, however, leads to Bloch oscillations which thwart the spec- 
tral flow. The momentum can escape only via relaxation processes. Taking these into account permits a 
physically transparent derivation of the mutual friction coefficients. 
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1) Introduction 

Imagine a two dimensional supcrfhiid, initially in its ground state, confined to the surface of a torus. 
Suppose now that a vortex-antivortex pair is created at some point on the surface and the vortex is moved 
slowly round one of the generators of the torus before being allowed to annihilate with its antivortex partner. 
One effect of this process is to give the superfluid order-parameter phase a unit winding number around 
the generator perpendicular to the motion of the vortex. The associated phase gradient implies that a 
supercurrent has been established in this direction. If no other momentum-carrying excitations were created 
along with the supercurrent, the system as a whole has acquired momentum perpendicular to the vortex 
motion. Moving the vortex therefore requires us to supply this momentum from an external source. This is 
the Magnus effect [1]. 

For a Bose superfluid this is all there is to the story: If we wish to move a vortex with respect to the 
background fluid we must (at least at low temperatures when there is no normal fluid component) place a 
wire or other object in the vortex core to supply the transverse momentum to the fluid. The reaction force 
the fluid exerts on the wire is the Magnus, or Kutta-Joukowski, lift force [2]. 

For a fermionic S'-wave superfluid* the situation is subtler because the vortex has low-energy bound 
states [3.4] whose role in the momentum balance equation has been studied for many years [5,6,7]. Recently 
Volovik [8] has cast a new light on this subject by showing that motion of the vortex with respect to the 
stationary condensate induces a spectral flow among these states. In a cartoon version of his theory this 
spectral flow generates, even with adiabatic motion of the vortex, a stream of unbound quasi-particles which 
carry off momentum equal and opposite to that of the induced superflow. The vortex can apparently be 
moved without any external source of transverse momentum. In this sense the spectral flow "cancels" the 
Magnus effect. The discrete nature of the bound state spectrum, however, complicates the picture. As 
observed in [8] and modeled in [9], a non-zero temperature is required to broaden the closely spaced levels 
so that they may behave as if the spectrum were continuous. It is only in the hydrodynamic limit that 
"cancellation" takes place. 

Despite the complicating necessity of level broadening, the spectral flow mechanism provides a very 
physical picture of the processes occuring in the vortex core. Several questions immediately arise, however. 

* We consider a neutral condensate for simplicity. The principal effect of the magnetic field in an Abrikosov 
vortex is to transfer the momentum supplied by the vortex to the positive ion lattice, thus ensuring that no 
superflow is induced beyond the penetration depth. Nothing significant changes in the core. 
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For example: What happens when a vortex is held stationary in a transport supercurrent? By galilean 
invariance, this situation is physically equivalent to a moving vortex and a stationary supcrfluid, yet - 
with no time dependence in the Bogoliubov-de Gennes equations — it is not immediately clear what drives 
the spectral flow. A second question is whether the spectral flow picture requires a modification of the 
conventional theories of momentum balance and Hall angle. The aim of the present paper is to discuss these 
issues within a simple model for the core states. 

We will introduce a quasi-classical picture, based on the geometric optics limit of Andreev scattering 
[10], for the evolution of the states in a stationary vortex core. This allows us to show that the discreteness of 
the spectrum leads to effects analogous to those in the one-dimensional Wannier-Stark ladder. Stark ladder 
resonances occur when a uniform electric field is applied to a Bloch electron in a periodic potential [11]. 
The electron initially accelerates but, in the absence of dissipation, is eventually Bragg reflected from the 
periodic lattice potential resulting in an oscillatory motion in a localized region. No net current flows - 
except that arising from the exponentially small inter-band Zener tunneling. If we describe this process in 
a gauge where A n = 0, A x = —Et we have a time-dependent hamiltonian and explicit spectral flow. In a 
gauge where A = —Ex, A x = there is no explicit time dependence but the same physics results. In both 
gauges dissipation via inelastic collisions allows the electron to avoid Bragg reflection and so permits a finite 
current. Analogously, relaxation processes in the vortex core allows some spectral evolution. By keeping 
track of the resultant momentum flux we find the consequences of the spectral flow for the vortex dynamics. 
These turn out to be the well-known mutual friction that couples the superflow to the normal flow via the 
vortex motion. The traditional Green function formalism of refs [6,7] must therefore tacitly take the spectral 
flow into account. 

The organization of this paper is as follows: First, in section 2), we exhibit a simple version of spectral 
flow and show how momentum entering the vortex core is recycled as quasi-particlc momentum. In section 
3) we will review the theory of the core states and its connection with Andreev reflection. In 4) we interpret 
the core state spectrum in terms of the failure of exact Andreev retro-reflection and, in 5), armed with the 
insight gained from this interpretation, we show how the spectral flow is mapped onto the Stark- Wannier 
problem. Finally, also in section 5), we account for the momentum flux to the normal component. 

2) One-dimensional Spectral Flow 

As a simple model of a vortex core consider the one-dimensional Dirac- Andreev eigenvalue problem 

-iv f d x A(x)e i9 ^ 
A(x)e- ie( ^ iv f d x 
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where A(x) = for < x < L (the "core") while A(x) = A = const, elsewhere. We take the phase of the 
order parameter to be 6{x) — 0l for x < and 6(x) = Or for x > L. Here vj is the Fermi velocity. We will 
use m to denote the fermion mass so that kf = mvf is the Fermi momentum and Ef = \vnv 2 is the Fermi 
energy. 



The bound-state solutions, & = 



, with e < A are easily found. The wavefunctions are of the form 



e + ikvf 
Ae~ i6R 

ae iex/v f 
fe-iex/vf 

e — ikvf 
Ae~ i9L 



e - k{x - L *> x > L, 

< x < L. 



(2.2) 



kx 



x < 0, 



with e 2 + (v/k) 2 — A 2 . Matching the solutions at x = 0, L fixes the ratio a/b and requires the eigenvalue e n 
to obey 

(2.3) 



For states deep in the gap, e <C A, this simplifies to 



(2.4) 



We see that if we gradually increase the phase difference across the core, A6 = Or — 6l, the entire 
spectrum moves up in energy. By the time AO has increased by 2tt each state has been replaced by the one 
below it. This is the spectral flow. 

The physical interpretation of the spectral flow depends on the context. If (2.1) were describing a charge 
density wave (CDW) system, the upper and lower components of W would be the amplitude of left- and 
right- going particles. Then a summation over occupied states gives the local charge density and current 

< &(x)V{x) >= < ip R {x)i) R {x) +i>l(x)i/) L (x) >=< p(x) > 

v f < * t (x)cr 3 *(a;) >=v f < {il) R {x)^ R {x) - i>{(x)i> L (x)) >=< j(x) > . 



(2.5) 



In a CDW a time rate of change of the phase of the order parameter induces a current < j >~ ^6, 
with the corrections being small when is small. Consequently the slow twisting of Or relative to 0l tells 
us that charge is flowing into the the gapless region < x < L. Since the time-dependent version of the 
Dirac-Andreev equation implies that p and j obey the conservation law 



dtp(x) + dxj(x) = 0, 
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(2.6) 



the inflowing charge must be accumulating in the gapless region [12]. Each time the relative twist increases 
by 2tt, a unit charge will have accumulated. In the same interval one of the (occupied) negative energy bound 
state levels has adiabatically crossed the zero energy level and taken the place of a positive energy state. 
The occupation number of the positive energy bound states has therefore increased by unity, consistent with 
the accumulation of unit charge. Eventually, with more twisting (the amount depending on L), the filled 
levels will reach the top of the gap and merge with the upper continuum. After this point each new unit of 
charge that flows in will appear as a low energy quasi-particlc. 

In a superconductor the upper and lower components of ^ are ip R , ip' L respectively. In this case the 
expressions for the current and charge density are interchanged 

v f < &(x)if/(x) >=v f < ^ R (x)tp R (x) +i> L (x)i/)l(x) > 

=vj < ip R (x)ip R (x) - ipl(x)tp L (x) >=< j(x) > 



= < ^ R {X)^ R {X) + ^ L {X)^ L { X ) > 
= < p(x) > . 



(2.7) 



The conservation law 



d t (ip R ipR -tpltp L ) + Vfd x (i/) R ip R + ip\ip L ) 

= a t * t * + v/a*o- 3 * (2.8) 

= o, 

also changes its physical interpretation as, on multiplication by kf, it becomes the equation of momentum 
conservation [13]. Now, instead of charge, each occupied bound state carries momentum The relative 

twisting of the phases on the two sides of the core represents an inflow of momentum from the condensate, 
and the spectral flow leads to its recycling as the momentum of low energy quasi-particles. 

We can make a simplistic model of the consequences of two-dimensional vortex motion by assuming that 
(to a first approximation) the process described in the introduction, the passage of a vortex around the L y 
generator of an L x x L y torus, can be mimicked as the breaking and reconnection after a 2tt phase twist of 
the order parameter in a collection of one-dimensional superfluids, one for each allowed k y = 2irn/L y . From 
the discussion following (2.8) we see that a single twist accumulates a momentum k x — ^J\kf\ 2 — k^ for each 
of the one-dimensional systems. This is the same amount of momentum we would get by translating the k x 
value of each particle on the torus by 5k x = ir/L x . If the total number of electrons is N, the net momentum 



accumulated in a passage of a vortex round the L y generator is therefore 

ttN irp 



m 



(2.9) 



where p = mN/(L x L y ) is the mass density. 

This is implies a rate of momentum accumulation in the core of 

dP x 
dt 



TT 

— pVy. 

m 



Since the circulation in a BCS vortex is n — Tr/m we see that momentum is accumulating, or being recycled, 
at a rate equal to the Magnus force on the vortex, F x — npv y . 

This cartoon version of the process is of course overly simplistic. The bound states in the two dimensional 
vortex core are not those of the one-dimensional equation (although we will soon see that they are closely 
related), and the variation of the order parameter in the y direction will couple the k y momenta so that they 
cannot be dealt with individually. In the next section we will begin to deal with these deficiencies. 



3) Two dimensional bound states 

In this section we will review the classical results of [3,4] on bound states in the core of a two-dimensional 
vortex. Our aim is to show that the physics of Andreev scattering reduces the full problem to a collection 
of one dimensional problems of the form considered in the previous section. 



We wish to solve the Bogoliubov-de Gennes equation 



3_V 2 
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(3.1) 



Here r and 6 are polar coordinates with origin at the vortex centre. 



For the moment we will leave the gap profile A(r) unspecified, but the angular dependence of the order 
parameter is such that the superflow is anti- clockwise with a single quantum of circulation. 

We now separate the radial and angular parts of the wavefunction. To do this we must first appreciate 
that u, v are invariant only under An rotations*. Therefore we seek solutions in the form 



u{r, ey e / 2 

v(r, 6)e' w ' 2 



ui{r)e ie / 2 



I e Z. 



(3.2) 



* The same is true for the one dimensional problem in section 2). There we saw that a 2ir twist in the 
order parameter shifts the particle momenta by 5k x = n/L x . If there were no quasi-particle created along 
with the twist this would lead to a double valued many-body wavefunction. Fortunately the x dependence 
of the quasi-particle serves to restore the single valuedness to the total wavefuction. 
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We find that ui(r), vi(r) obey 
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(3.3) 



To separate the rapidly varying degrees of freedom from the slow, we follow [3] and write 

H£\k f r)+c.c 







fl 


Vl _ 




.91. 



(3.4) 



where Hm(kfr) with m = \Jl 2 + \ is a Hankcl function, and /j, gi are slowly varying envelope functions. 
We may approximate the Hankcl function by the WKB form 



(3.5) 



which is valid for r larger than the WKB turning point r = b. This turning point is related to kf and I by 



Substituting (3.5) in (3.3) and keeping only the large terms we find 

-ivf Vr ' 2 - b ' 2 d r 

A ivf- 



A 

■ Vr ' 2 - b ' 2 d r 





'fl' 




'fl' 
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.91. 



(3.6) 



We can make (3.6) more intelligible by trading r for a variable x defined by r 2 = x 2 + b 2 [4]. We find 

I 



-ivfd x 
A 



A 




7" 




iv f d x _ 
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2m(x 2 + 6 2 ), 
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7" 




.5. 



(3.7) 



This substitution has a geometric meaning: For any Bessel function Ji(kr) the turning point distance, b, 
is the impact parameter of a particle moving past the origin in a straight line with momentum k and angular 
momentum I. The quantity x is the distance of the particle from the point of closest approach to the origin. 
Along with x we introduce a polar angle "0" by x — btan9. This angle differs only by a constant from the 
true polar angle 8, so we temporarily abuse notation and make no distinction between the new angle and 
the original. We then notice that 



I 



_ _ 1 <W_ 

2m{x 2 + b 2 ) ~ 2 Vf dx 



so that, defining, 



'fl' 




\ fie l0/2 1 


.91. 




gi^l 2 _ 



(3.8) 



(3.9) 



to undo the transformation that removed the angle dependence from the order parameter in (3.1), we get 

T f fo Ae !l =*V-]- (3-10) 
Ae w iv f d x \ [gi\ [gi J 



This looks very much like the one dimensional eigenvalue problem solved in section 2). It is not yet 
identical, however. In (3.10) the coordinate x — \Jr 2 — b 2 is restricted to positive values. Furthermore 
boundary conditions have to be imposed on fi , gi at x = to ensure that the Hankel functions can combine 
to give the J l±k (jfer) = \{H^ L {kr) + H^{kr)) which are finite at the origin. We may nonetheless extend 

2 ^ t± 2 t± 2 

x to negative values by regarding the part of (3.4) with the incoming Hankel function Hm (kjr) as living 
on the negative x axis, and the outgoing part Hm\kf-r) on the positive x axis. With this interpretation the 
boundary conditions at x = translate into the requirement of continuity of fi, gi there. 

Physically the transformation of the two-dimensional eigenvalue problem (3.1) into the one-dimensional 
(3.10) occurs because each bound quasi-particle is bouncing back and forth along a straight line, its direction 
of motion being repeatedly reversed by Andreev scattering off the increasing value of the gap. 




Fig. 1. A bound state with an electron (solid arrow) being Andreev reflected as a hole (dashed 
arrow) . 

The A(r) profile found from a self-consistent solution of the Bogoliubov-de Gennes and gap equations 
has scale R sa Vf/A [14]. In order to have analytic expressions for the bound state eigenvalues, we will not 
use such a self consistent A(r), but take instead a step function: A(r) = 0, r < R and A(r) = A for r > R. 
We will also assume R to be somewhat larger than Vf/A so that we can ignore, in (3.10), the variation of 
9 in the regions where the wavefunction is evanescent. With these assumptions we can directly apply (2.3), 
(2.4). 
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Each quasi-particle trajectory now coincides with a chord of a circle of radius R. If this chord subtends 
an angle of 2% at the centre of the circle it has length L = 2i?sinx (see Fig. 1.). Furthermore the quantity 
A9 = Ol — Or, the difference in order parameter phase at the two ends, is equal to 2%. For states with energy 
e <C A the energy eigenvalues are therefore 

en(0 = 7p^ (2x + 27r(n+l)). (3.11) 
AR sin x 2 

Here I — bkj = Rkj cosx- For small / and n = — 1, (3.11) reduces to 

e_i(0 = -w O i, w "=2^- (3J2) 

This is the topological branch of low energy excitations which is important for the spectral flow. Note that 
u) = (2mR 2 )~ 1 is the angular velocity of the superfiow at the boundary of the core. 

The maximum value of I occurs where the chord length becomes zero and e — ► A. This maximum value 
is ^max = Rkf ~ Vfkf/A = 2Ef/A w 10 3 — 10 4 . This is large enough that I can almost be regarded as a 
continuous parameter. We will often use this observation to write expressions such as de/dl without further 
comment. 

In the next section we will show that the I dependence in (3.11), (3.12) is a consequence of the failure 
of the Andreev scattering to be perfectly retro- reflective [15]. 

4) Andreev Reflection and the Bound State Spectrum 

The bound states may be thought of as standing waves set up by Andreev reflected quasi-particles 
repeatedly traversing a chord of a circle. For a bound state of definite angular momentum, I, the orientation 
of this chord will be indefinite. To produce a state localised on a chord with specified orientation and impact 
parameter we must take a linear combination of angular momentum states, 



>= 



5>e- iW |i>, (4.1) 



where the coeficients a; are large only for I w b/kf. Because of the Rayleigh criterion relating the extent of 
a plane wave-front and its diffraction spread, there will be an uncertainty relation between the sharpness of 
6 and b. 

We may follow the time evolution of these wave packets by using the small I approximation e_i(Z) = 
— lo^I. We find 

\6,t >=Y J a l e ~ iWe+iluJot \ l >= \{0-oj o t),t = > . (4.2) 



The chord is thus seen to precess at the angular frequency ujq in the sense opposite to the superflow. If we 
use the more general form (3.11) for the energies, we will need to replace — ujq by a group angular velocity 

de(l) 



u> 



dl 



(4.3) 



but the same general picture will hold. In this section we will sec that this precession can be understood 
rather precisely as the combined effect of two processes each of which causes the Andreev reflection to fail 
to be be perfectly retro-reflective. One of these processes is due to the superflow and we will be able to 
use this insight to determine the evolution of the bound states when the vortex is immersed in a transport 
supercurrent. 



To begin with we consider a plane SNS junction of width L having a supercurrent v s flowing parallel to 
the junction (see Fig. 2.). 



V. 




Fig. 2. Bound state in an SNS Junction with superflow. The solid arrow is the electron trajectory 
and the dashed arrow the hole. 

There will be bound states in the junction with momenta close to k = {k x ,k y ) — (kf cos#, kf sin( 
These states arc found from the Andreev hamiltonian 



H = 



-i(k x d x + k y d y ) A(x)e 2imv ° y 

A{x)e~ 2 " nV °y i(k X 8 X + kydy) 



where A(x) = for < x < L and A(x) = A elsewhere. They have the form 



i(Sk^x+8k y y) 
e — i(Sk x x+8k y y) 



< x < L. 



(4.4) 



(4.5) 



Here Sk y = mv s and (for states with e <C A) Sk x — jr[2ir(n + |)]. Their energy is 



e = Vf(Sk x cos 8 + Sk y sir 



(4.6) 
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In the figure we have drawn the trajectories of the electron and hole as if they were at a definite location 
in the junction. In reality, specifying a y location requires making a wave packet 



\y>= J ^-a(k)e- ik *y\k>, 



with some spread of k y . This packet will drift in the +y direction at the group velocity 

de(k) tan^ 



Wdrift = 



<9fc„ 



m 



+ v s 



(4.7) 



(4.8) 



For trajectories that are at close to normal incidence on the supercurrents (9 = or k y w 0) this velocity is 
essentially v s . 




k+5k 



Fig. 3. The electron (solid arrow) and hole (dashed arrow) momenta do not lie exactly at k on the 
Fermi surface, consequently their momenta differ in direction by a small angle 59. 



This drift can be accounted for by noting that the Andreev scattering is not perfectly retro-reflective. 
Consider reflection from the right-hand supercurrent in Fig. 2. Although both the incident electron and the 
reflected hole have nominal momentum k, they in fact have momenta k± 5k. (see Fig. 3.). 



This means that each reflection causes a small change in the angle of incidence 

59 =2(-5k x sin 9 + 5k y cos 9) 

= _ tan0+— i -. 

Ef Vf cos 6 



(4.9) 



Because both particle and hole move with speed Vf, it is easy to see (see Fig. 3) that such a change in angle 
leads to the trajectory migrating up the junction with velocity 



v f 

Udrift = J O. 

2 cost' 
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(4.10) 



Inserting (4.9) in (4.10) precisely reproduces (4.8). 




Fig. 4. The slight difference between the angle of incidence of the electron (solid arrow) and the 
angle of reflection of the hole (dashed arrow) causes the bound state to migrate up the SNS junction. 



The same result holds true for the bound states in the circular core. From section 3) we know that the 
bound state energies are 

(4.11) 



where I = kfRcosx- The group angular velocity is therefore 



de cos x e 



dl sin %fc/J? 2mR 2 sm x 



(4.12) 



For states with small e (i.e. n = —1, x ~ 7r /2) the latter term dominates. 

The geometric optics picture works here also. We can use the plane interface formula for 89 after noting 
that the angle of incidence, called 9 in (4.9) is now 6 = \ — X- Notice, though, that a positive 89 means that 
the point of impact of the return trajectory moves backwards through an angle of 289 (See Fig. 5.). This is 
the origin of the minus sign in e(l) w — ujqI. The other infomation we need is the time between impacts at 
the points A and C. This is 8t = ARs'mx/vf. Putting all the ingredients together gives 



cosx e 



1 



sin 2 x^/-R 2mi? 2 sin 2 



(4.13) 



as before. 
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Fig. 5. Andreev reflection from the supercurrent causes the chord defining the bound state to 
precess in the opposite sense to the superflow. 

The I dependence of our core-state spectrum is therefore entirely determined by the geometry of Andreev 
reflection. The only role of the wave character of the states is in the condition that I be an integer, and in 
the n dependence of the spectrum. Since only the zero-crossing branch, n = — 1, is of interest in the spectral 
flow problem, the n dependence is irrelevent here. 

The lesson we learn from this section is that the quasi-classical motion of the chord is governed by two 
processes whose effects add. One, proportional to the bound-state energy is independent of the superflow, 
and the other, proportional to the superflow, is independent of the energy. The additivity of the effects also 
applies when the superflow is not parallel to the SN interface. We will need this fact in the next section. 
5) Frustrated Flow and Mutual Friction 

If we insert our vortex in a transport supercurrent by replacing the Ae l6 gap function by Ae l6 e 2lmVsX 
then either considerations of the group velocity or of the geometry of Andreev reflection will show that the 
wavepackets will try to migrate across the core at speed v s , just as in the parallel sided junction. In the 
absence of precession, this motion across the core will cause a monotonic change of [6r — Ol) from (where 
e = —A) to 2n (where e = +A) leading to a steady stream of states crossing the gap and becoming unbound 
quasi-particles, just as in section 2). 

This migration is, however, in competition with the — u>o precession. The competition can be described 
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quantitatively by writing an effective hamiltonian governing the evolution of the wavepackets 

H = -ujqI + v s • k, (5.1) 

where k = (fc/ cos#, kf sin#). The observables I and 8 are canonically conjugate and have commutation 
relations [8, 1] = i. 

For a supercurrent in the +x direction, the orientation and impact parameter of a wave-packet of bound 
states evolve according to the quasi-classical equations 

• oh 

dH (5 ' 2) 
/ = - — =k f v s sin9 

The orbits are 

I = cos wot + const. (5.3) 

We see that a packet with an initially positive value of sin# tries to cross the core from — |£ max | to + |? max | 
but, because of the evolution of 0, it is reflected back. This is not surprising since the hamiltonian (5.1) is 
identical to the Wannier-Stark hamiltonian for a tight-binding model. In this interpretation I would label 
the atomic site, wo would be the external electric field and 8 the crystal momentum. The 2tt periodicity of 
the latter is a direct consequence of the discreteness of the values of I — just as in the vortex the discreteness 
of I is a direct consequence of the 2ir periodicity of 8. The to-and-from motions of l(t) are therefore Bloch 
oscillations caused by the discreteness of the bound state spectrum. 

The only way a state can make the journey from — |Z max | to + |^ max | without being reflected would be 
for v s kf /u>q to be larger than Z max . But this is equivalent to the condition that v s exceed the pair-breaking 
velocity, v pa-il k f = A. 

Although the spectral flow is thwarted by the Bloch oscillations there are still effects from the evolution 
of the states. We can analyze these by writing a Boltzmann equation for the distribution function, n(l,8), 
in the 1,6 phase space. We have 

dn dn . dn 

--^- + kf v s s m 8- = I (5.4) 
where / is the collision integral. We will use a simple relaxation time approximation to (5.4) by setting 

I=--(n(l,0)-n(l,6)) (5.5) 

T 

where fi(l,0) is an equilibrium distribution whose exact form is determined by the normal velocity v„. (In 
a superconductor the normal component relaxes to the lattice, so in this case v n = lattice-) For example if 
v n = v s = we would expect n — > no = 8(1) so that all negative energy states are occupied. 
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It is actually more convenient to write equations for the moments of n that give the total momentum 
in the core 

< k x >=- ( dl — (n(l,9) - n (l,6))kf cos 6 

2 J 2?r m a\ 

1 f d9 ( ' 

< k y > = - / dl — (n(l,6) - n (l,6))kf sin 6. 

2 J 2n 

The factor of \ in front of these integrals compensates for the double counting of particles and holes, and 
no(l, 9) = 9(1) reflects the normal-ordering which ensures that < k >= when everything is at rest. 

As an example suppose we fill all the negative energy states of (5.1) by setting n(l, 9) — 9(oj l~v s kf cos 9). 
We find < k y >= and 

1 k2 f Vs v s , „. 

< k x >=---!— = -pn— 5.7 
4 ua lu 

where n = ir/m is the quantum of circulation and p = mk'j /Air is the equilibrium mass density of the electron 
fluid. This would be the appropriate equilibrium distribution for the case v„ = 0. The non-zero value of 
< k x > reflects the reduction in current through the vortex core because of the non-zero value of p n there. If 
v„ = v s , on the other hand, we expect < k >= because, even though p n ^ in the core, the total current 
j = p s v s + p n v n is unaffected by the vanishing gap. 

Using these definitions and insights we find 



d<k x > 1 

— =w < k y > +npv sy - - (< k x > - < k x >oj 

d < k v > 1 

loq < k x > -npv sx (< k y > - < k y > ) 



(5.8) 



dt t 
where 

<k> =^(v„-v s ). (5.9) 

Wo 

In the steady state flow < k x . y > will be constant, so the left hand side of (5.8) will be zero. We can 
therefore solve for < k XtV > and find the rate at which momentum is being transferred to the lattice 

" k ' M ' : k • k •, i. (5.10) 



dt 



lattice 



T 



We find 

' d < k x > \ loqt 1 



— I — Kp- — ■ — j V nx + Kp— — ^V ny 

dt /lattice l+^o T 1+^T 2 

d < k v >\ ujq t 1 



dt Aattice" ^l+^y ^l+CoV 

Notice that v s does not occur in (5.11). 
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These calculations have been made with the vortex stationary. Because of the galilean invariance of the 
overall system, we may replace v„ by (v„ -vj,) when the vortex is moving at velocity v^. The rate of loss 
of momentum from the core to the lattice is then 



d < k > 
di 



= D{w L - v n ) + D'i x (v L - v„) (5.12) 

lattice 



The quantities D, D' are the mutual friction coefficients. Their values 

LOqT 



D up — — 
1 + ujAt 2 



D' =up 



n ' 

(5.13) 



1 + wgr 2 



are the same as found, for example, in [16]. For the purposes of making comparisons with other work we must 
point out that in deriving these expressions we have assumed that we are at sufficiently low temperatures that 
there are very few quasi-particles outside the vortex core. We have therefore made no distiction between 
p and p s . Similarly the linearization inherent in our use of the Andreev equations makes it difficult to 
distinguish between p and the constant Co which is defined as the electron mass times number of states 
within the Fermi sphere, and so coincides with the normal-state density. 

Equation (5.13) determines that Hall angle. While the task of finding a solution to the full dynamical 
gap equations for a vortex immersed in a superflow is rather daunting, we do know that such a solution must 
satisfy the law of momentum conservation. The vortex core velocity Vl must therefore satisfy the momentum 
balance equation 

= npz x (v L - v s ) - D(v L - v„) - D'z x (v L - v„) (5.14) 

as a consistency condition. The first term is the Magnus force which gives the rate at which momentum 
enters the core. The terms with D and D' give the rate at which momentum is being lost from the core to 
the lattice by mutual friction. 

There are two extreme cases. In the collisionless regime, wt > 1, the vortex has no choice to move 
with the superflow. In the opposite, hydrodynamic, regime cj t is small. Then the Vl part of the D' term 
almost cancels the Vl part of the Magnus force term. This allows the vortex to move at right-angles to the 
superflow, unwinding the order parameter phase gradient and dissipating the superflow. 

It is worth pointing out that the cancellation between the two vl terms in the hydrodynamic limit 
should not be regarded as a "cancellation" of the Magnus force. To claim that it is, is akin to saying that 
the the "F" cancels the "ma" in F = ma. This may seem a mere quible, but an an inappropriate choice 
of language has caused confusion in the literature. It is also worth stressing that while the reactive term 

16 



D'z x (vf, — v„) does not cause any dissipation of energy, its presence is entirely due to quantum inchoherent, 
entropy-generating relaxation processes. 

6) Conclusion and Discussion 

We have seen that a quasi-classical geometric optics model gives a good description of the processes in 
the core of a vortex immersed in a transport current. Even though there is no explicit time dependence in 
the Bogoliubov-de Gennes equations, spectral flow would still occur were it not suppressed by an analogue of 
Bloch oscillations originating in the discrete nature of the spectrum. The spectral evolution does, however, 
lead to a non-equilibrium occupation of momentum-carrying core states. When we account for the processes 
by which the occupation distribution attempts to relax, we find expressions for the mutual friction parameters 
which coincide with those obtained by traditional Green function methods [6,7,15]. It is therefore clear that 
the Green functions tacitly include the physics of spectral flow. The new method of explicitly following the 
evolution of the states does however have the advantage of being much more physical. 

Our general conclusions are consistent with recent work by Kopnin et al. which is briefly reviewed 
in [17]. Some other recent discussions of the Magnus force and momentum balance are by Hoffmann and 
Kummel [18], Gaitan [19], Ao [20], and Simanek [21]. The first of these papers has much in common with 
the present approach. In particular these authors focus on the role Andreev scattering plays in transferring 
momentum from the condensate to the core states. The next two use the Berry phase approach to confirm 
that the Magnus force correctly gives the momentum flow into the vortex core. The last combines effects of 
external and core states to find an effective "Magnus" force that is a combination of the true Magnus force 
and the mutual friction terms. 
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